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0^ ■ ABSTRACT 

I ■ A convenient form of kinetic equation is derived for pair annihilation of heavy 

D I stable particles relevant to the dark matter problem in cosmology. The kinetic equa- 



tion thus derived extends the on-shell Boltzmann equation in a most straightforward 
way, including the off-shell effect. A detailed balance equation for the equilibrium 



> 

X. 

H , abundance is further analyzed. Perturbative analysis of this equation supports a 

previous result for the equilibrium abundance using the thermal field theory, and 
gives the temperature power dependence of equilibrium value at low temperatures. 
Estimate of the relic abundance is possible using this new equilibrium abundance in 
the sudden freeze-out approximation. 



1 Introduction 

Generalization of the Boltzmann equation is important both from the point of 
fundamental physics and of application to many areas of physics. Since its investi- 
gation involves quantum dynamics at finite time, its general formulation in a form 
of practical use is undoubtedly difficult. One important classs of physical situations 
however seems amenable to relatively straightforward analysis. It is a small system 
immersed in a large thermal environment. 

Application of this class of physical situations includes pair-annihilation of heavy 
stable particles in cosmology; the dark matter problem [|l| . The conventional practice 
for estimate of the relic abundance of dark matter particles is first to compute the 
freeze-out temperature below which the annihilation effectively ceases due to rapid 
cosmological expansion, and then to use the ideal gas form of number density at 
that temperature. A thermally averaged Boltzmann equation is a basis for such 
computations . Clearly, the use of the on-shell Boltzmann equation and the ideal 
gas distribution function should be questioned, especially at low temperatures such 
as T < M/30 (M being the mass of the heavy particle) with very small e~^'^l'^ like 
10^^'^, which happens to be the relevant temperature region for the dark matter 
annihilation. 

As in our previous studies [@],@],[|1, we take in the present work of quantum 
kinetic equation the method of integrated-out environment in much the same way 
as some treatment of quantum Brownian motion |^. This way one can well cope 
with quantum mechanical behavior at finite times beyond the infinite time behavior 
taken care of in the S-matrix approach of the Boltzmann equation. In this paper we 
report a great and practical simplification over our previous published result [^. We 
are thus ready for a realistic study of supersymmetric dark matter problem. 

The equilibrium number density thus computed agrees with that obtained using 
the thermal field theory in . The result of the present work is however more general; 
we present the full quantum kinetic equation, not just the equilibrium result. The 
essential ingredient in this formulation is the Hartree approximation, which makes 
it possible to use the result of the solvable model [Q, with an extra assumption of 
slow variation of physical quantities. 



2 Model of pair annihilation 

Suppose that a pair of heavy particles ip (boson in this case) can annihilate into 
a X (boson) pair; ifip -^ XX with a dimensionless interaction strength A. A discrete 
symmetry under ip -^ — y? is imposed to forbid the p decay, p -^ xX- We assume 
that the lighter particle x makes up a thermal environment of temperature T = 1/(3 
in our unit of the Boltzmann constant ks = 1 ■ Thermalization becomes possible 
either due to their own self-interaction characterized by a dimensionless coupling 
A^, or interaction with other light particles which need not to be specified for our 
purpose. A relativistic field theory model we consider is thus given by the Lagrangian 
density 
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C^ + C^ = -{d.^r - -M' p' + -{d,xr - 7;m'x' , (2) 

C-... = -^v'x'-^v'-^x' + SC. (3) 

The last two terms (oc A^ and A^) in the interaction are introduced both for consis- 
tency of renormalization and for thermalization of x particles, hence the parameters 
are taken to satisfy |A^| ^ A^ ^ 1 , |Ajj,| < 1 , but |A^| ^ |A|. The renormalization 
counter term 6C is given by 

'-^is.'Pr - ^/ + '-^(a^xr - ^x= - '-^^V - ^/ - '-^x' . (4) 

We focuss on a particular dynamical degree of freedom, the heavy p field, and 
integrates out the environment part, the x field, altogether. In the influence func- 
tional method [§] this integration is carried out for the squared amphtude, namely 
for the probability function, by using the path integral technique. This way one has 
to deal with the conjugate field variable p' along with if, since the complex conju- 
gated quantity is multiplied in the probability. The influence functional JF is thus 
defined by 

J '^^ J ^^' ^^P ^ J ^^ i^xi^) - ^x'i^) 

+ Cinti^{x),xix))- Cintiv'{x),x'ix)))\ . (5) 

We convolute with the influence functional the initial and the final density matrix 
of the X system. For the initial state it is assumed that the entire system is described 



by an uncorrelated product of the system and the environment density matrix, 

p, = p('^) X pS^) , pI^) = pf^ = e-/5^«(^)/tre-^^°(^) . (6) 

Here p^- is the density matrix for a thermal environment. The density matrix pi' 
for the (fi system is arbitrary. The choice of the initial density matrix pi is not crucial 
and the final result is expected to be insensitive to this choice provided that it does 
not commute with the total hamiltonian, [pi ,H]^0. We need this condition to 
guarantee a departure from complete thermal equilibrium. At a final time tf the x 
integration is performed taking the condition of non-observation for the environment, 



dXf J dx'f6{xf-x'f) {■■■), (7) 

with the understanding that the environment is totally unspecified at the time tf. 

The result of x integration is given by a series of Gaussian integral if we expand 
in powers of A^ of the A^ x^ interaction, and the influence functional to order A^ is 
of the form, 

exp[-- / dxdy {^2ix) a R{x,y)^2iy) + i ^2ix) a i{x,y)X2{y))], (8) 
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X2{x) = ip'ix) + ^'\x) , U^) = ^'{x) - ^' \x) , (9) 

a{x , y) = anix , y) + iai{x , y) 
= A^ (tr {T[x\x)x\y) P^^'^]) - (tr xVj.^))^) , (10) 

if one replaces the mass term in the original Lagrangian by a 0[A] temperature 
dependent mass, M'^{T) = M^ + ^ ^^ valid for m <^T. Note the presence of the 
time ordering, xq > yo, in the above formula. The kernel function a{x,y) satisfies 
the time translation invariance, thus may be written as a{x — y). An explicit form 
of the kernel function a{x) or its Fourier transform is given later in eqs. (p3D and 
(|25|). Higher order terms in A^ are actually present in the exponent of the influence 
functional. These contribute either to many-body processes we are not interested in, 
or to negligible higher order terms to our process. 

The convolution with the heavy ip system gives the reduced density matrix at 
time tf] 

(11) 



from which one can deduce physical quantities for the ip system. Here S{ip) is the 
action for the if system obtained from the basic Lagrangian. 

The correlator is defined using the reduced density matrix; for xq > Vq 

{ip{x)ip{y)) = J dip{x) J dif'ix) J d(p{y) J dip'{y) 5 {ip{x) ~ (p'{x)) 

V^ J V^' e^^(^)-^(^') J-4[v. , y.'] ^{xMy) p'^''\v{y) , ^\y)) . (12) 

The path to be taken here is those defined on the time interval, xq > t > y^. 
A nice feature of this correlator formula is that the initial memory effect appears 
compactly via the reduced density matrix p^^\ With the Heisenberg evolution 
^{x) = e'^^'^ip{x ,0)e~^^'^, this correlator is equal to 

ti Pi ip{x)ip{y), (13) 

where pi is the total initial density matrix (^. Whenever this density matrix does 
not commute with the total hamiltonian, [pi ,H]^0, this correlator does not respect 
the time-translation invariance, 

{^{x,ti)^{y, ta)) ^ {^{x,ti - t2My, 0)) , (14) 

unlike the correlator in complete thermal equilibrium. 

The model thus specified is difficult to solve due to the appearance of the quartic 
term oiip in the influence functional (^. The situation is however simplified when one 
considers a mean field approximation. In the mean field or the Hartree approximation 
one replaces an even number of field operators by a pair of two operators times the 
rest of several averaged two-body correlators. This approximation is good if one 
can ignore a higher order correlation than that of two-body. The Hartree model is 
expected to work well in a dilute system characterized by a low occupation number 
for each mode. 

The Hartree approximation we introduce is then a Gaussian truncation to the 
infiuence functional; we replace the original one by properly defining new real kernel 
functions Pi{x ,y) in the quadratic form; 

^2y,^']=exp[- / dxdy {^{x)PR{x,y)^{y) + i^{x)Pi{x,y)X{y))], (15) 
X{x) = <f{x) + <f'{x), ax)=fix)-v'ix). (16) 

When the time-translation invariance holds for /5j, one can work out its consequences 
to all order of A for a given kernel (3i{x,y), since it becomes a solvable Gaussian 
model. 



For derivation of the self-consistency relation introduce the full propagator by 

G{x,y) = i{^{x)<^{y)):F=:Fi. (17) 

using the influence functional (^. We compute correlators, 

(X(x)e(y)) , {X{x)X{y)) , i^xny)) , (18) 

in two theories, using the two different forms of the influence functional, JF4 and the 
truncated one J^2, and identify two results. One then has {P = Pr + iPi)] 

j3(xi , X2) = —i a(xi — X2) G{xi , X2) ■ (19) 

It is important to note again that there is no translational invariance with respect 
to time variables ti,t2, since the system is not in complete thermal equilibrium. 
Thus, one should distiguish two time variables for the two-point correlator, the 
relative time t = ti — t2 and the the central time t = [ti + t2)/2. A fundamental 
strategy in the present work is to note the slow change of physical quantities for 
the central time and exploit it in the analysis of the above relation (0). One thus 
solves short time dynamics of the relative time assuming a constant central time. 
The result of solvable Gaussian model 0], [Q may then be used provided that the 
kernel function f3 is given. Since the correlator thus derived is a functional of an yet 
unknown /?, one arrives at the self-consistency equation for this (3, which contains 
the constant central time. This constant central time is then allowed to slowly vary 
in the spirit of adiabatic approximation. We shall explain this procedure in detail in 
due course. 



3 Correlation function in Hartree approximation 

It is convenient to define the Fourier transform with respect to the relative vari- 
able, 

a{k, t)= I d\x - y) (^(x) v^(y))e^^-(^-^) . (20) 



Since spatial homogeneity holds, this quantity a may depend on the central coordi- 
nate only via t = (xq + |/o)/2. 



Suppose that the central time dependence is slow, and we may take this as a 
constant. The correlation function can then be calculated using the technique of 
the generating functional for the influence functional, as explained in ref. ||^. In the 
Gaussian model or in the Hartree model of pair annihilation each Fourier component 
of correlators is given in terms of a basic function g{k , t) and the density matrix of 
the initial state. Suppressing all Fourier momenta and denoting q{t) for Fourier 
component of field and p{t) for its time derivatives, one has, for instance, for the 
q — q correlator 

{q{ti)q{t2)) = - J gih - h) + r dt r ds g{h - t)f3it{t - s)g{t2 - s) 
2 JO Jo 

+ 9{ti)g{t2) pi + g{ti)g{t2) g? + ^ ( ^(^1)^(^2) + g{ti)g{t2) ) Tp^^ ■ (21) 
The initial density matrix elements are needed only to sepcify three expectation 



values, p1 , qf and {pi , g* }. 

The basic function g{t) is related to the spectral weight r{uj), and ultimately to 
the Fourier transform of the correlator. Fourier transforms, r and r^, with respect 
to the relative coordinate are introduced as 

/3(x,l/) = /^r(fc,:^)e-M.-), (22) 

One obtains from the self-consistency relation 

2 [^ d^k' r^{k-k') 



r(k, t) = 2y / — — - "^ ^,^ \.,M k\ t) . (24) 



The important relation between r{k ,t) and a{k,t), eg. (|2^ , is a result of eg. (p!9|) . 
The 3-momentum dependence is always even, being function oi \k\ = k. Thus, r^ is 
a function of two variables, k^ = uj and fc, hence can be written as r^(u; , fc). It is 
calculable from Fig.l Q] and is odd in uo; for cj > 



rJuj , k) = —^ 1/ 1 - -4^ diuj - Vk^ + 4m2) + ^ In - — ^---- , (25) 
^^ ' ^ 167r2 I V cj2-A;2 ^ ^ ^ /Sk 1 - e-^l-^-l ' ' ^ ^ 



^± = 77±oJl--^^- (26) 



u! k / 4m2 

2 2 V ^ ~ ^^^^^ 
The basic function g that appears in the correlator is then given by 

(■00 _ 



gt{k,T)=i duH{uj,k,t)e-"^\ (27) 



00 



^^^'^^ (P-M2(T,t)-n(fc,t))2 + (7rr_(fc,t))2' ^^^^ 

Uik,t)^vr^''-^^, (29) 
where even and odd parts of the spectral weight are defined by 

r±{k,t) =r{k,t)±r{-k,t). (30) 

The time dependence of the mass shift M'^{T , t) is weak and its asymptotic equihb- 
rium value M'^{T) to 0[A^] is given later by eg . (|55|) . With the help of 

ht{u;,k,T)= r dsg,{k,s)e-''-^= T rf^^ ^^"'^' ^ '/^ (e^'^ - e^) , (31) 

JO J-oo LU — U 

the correlation function is calculated as 



{(f{x)if{y)) 
d^k 



yt{k,x'-y')+ r dco'-±^^^^h;ico,k,x')h,ico,k,y') 

Z J —OQ Z 



(27r)3 
+ gt{k, x^)gt{k, y°)^ + gt{k, x'>)gt{k, y°)g^ 

+ ^ {9t{k, x^)gt{k, /) + gt{k, x^)gt{k, V^)} <ik(l-k + W-k 



^.k-i^x-y) _ ^32) 



Note that t = (xq + yo)/2. The integrand of 3-momentum integral of this formula 
is divided into the antisymmetric part (1st term) and the symmetric part (the rest) 
under the interchange Xq <-* yo- On the other hand, the antisymmetic part of the 
correlation function is related to the spectral weight H as seen from eg. (pT]) , thus 
the antisymmetric part of the correlator is given by 

a.{k,t) = 27rH{k,t). (33) 

The eguation thus derived (^) can be regarded as a self-consistency eguation 
for {<f{x) <f{y)), or eguivalently for its Fourier transform a{k , t). Its derivation rests 
with the Hartree self-consistency and the slow variation approximation on the cen- 
tral time. In the following sections we derive more convenient form of differential 
eguations that gives an equivalent result. 



4 Detailed balance and equilibrium abundance 

One takes the infinite central time limit of the self-consistency equation in the 
preceeding section to derive detailed balance relation, from which one can analyze 
the equilibrium abundance. 

First, one has for the correlator 

in the infinite central time limit, (xq + yo)/2 —>■ C)0. All initial memory effects drop 
out in this limit, as expected. Its Fourier transform gives odd and even parts, 

a_(fc,oo)=27ri7(A;,oo), a+{k, oo) = ^^^+(^' ^^ H{k, oo) , (35) 

r_[k, cx)) 

since the spectral function H{k, oo) is odd in ko- 
We further define a quantity, 

rik,^) = '-t^. (36) 

cr-{k, ooj 



One set of the consistency equation derived by combining eqs. ( PSJ ) and ( P^ ) 

r{k, oo)T{k, oo) — r{—k, oo)(l + r(A;, oo)) = , (37) 

is then written as 

2 f^ dpQ f d^p' f d^ki f d^k2 ,3 



{ 5(4)(p + p'-k,- k2) K(l + /i)(l + h) - (1 + r)(l + r')/i/2] 
+26^'\p + p +k,- k,) [rr7i(l + /a) - (1 + r)(l + r')(l + /i)/2] 
+5(4) (p + p' + fc, + A;^) [rr'hh - (1 + r)(l + r')(l + /i)(l + /s)] 
+5(4) (p _ p' _ A;i - A;^) [r(l + r')(l + /i)(l + /2) - (1 + r)r7i/2] 
+2^(4) (p ^p' + k,- k^) [r(l + r')/i(l + /2) - (1 + r)r'(l + h)h] 
+5(4) (p _ p' + A;i + A;^) [r(l + rO/1/2 - (1 + r)r'(l + h){l + /a)] } . (38) 

Here r = r(p , 00) , r' = t{jp' , 00) , while 



Uk = ^k'^ + m\ (39) 

fi = U{k,), f2 = U{k2), U{k) = ^^J _ ^ . (40) 



Note that the momentum p is not integrated in this equation, thus eq.([38|) is an 
integral equation for the unknown function r(p , oo). 
There is an obvious solution for this set of equations: 



due to the presence of energy conservation in each associated process of eq. (|38|) . The 
3-momentum dependence in r(p , oo) is missing. We further advance our analysis 
assuming that this is the unique set of solutions to the above equation. Putting 
these into the other set of consistency equation at infinite time, one obtains a closed 
form of self-consistency equation for r_; 



/oo W^t-' 
-— ^ i7(A;' , oo) r^(fc + fc' 



(27r)3 ^'"'^^>^^ ^(e/3(fco+fc'o) _i)(i_e-/3fc'o) ' 

(42) 

^^^'°°'' " (P - M\T) - n(A;, oo))2 + (7rr_(A;, oo))^ ' ^^^^ 

The quantity J\{k , oo) is given by a principal value integral of r_(A; , oo) like (^9]). 
The other independent quantity r+ is determined by 

r+(fc, oo) = coth r_(A;,oo), (44) 



hence giving the correlator 



d^k 



This way one can completely determine r(A; , oo), hence cr(A; , oo). 

One may first note the trivial case of application of eq. (^Sf) , taking the zero 
coupling limit r = 0. In this case one should understand the quantity H as 

H{k,oo)=6{kl-El)e{ko), (46) 

with 



Ek = yk^ + M^, (47) 

thus gets the free field correlator at finite temperature T. 

Next, perturbative analysis of the detailed balance relation shows that to 0[A^] 

d^k' 



r_(fc,oo) = 2A2 f 



{27ifEk 
+ r^{kQ - Ek,k- k' 



r^{ko + Ek,k + k' 



g/3A:o 



(^^fi{ko+Ey) _ l)(l_e-/3^fc') 
g/3fco _ 1 
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(48) 



The spectral function r_ given by ( ^S]) may be written in a suggestive form by using 
an expression for r^ of 0. It is equivalent to 

r_(p,oo) = — f dUp, f dUk f dnk'i27rY5^^\p + p' + k + k') 

{ 5{po + Ep, -ujk- ujk') [/p'(l + /fe)(l + fk') - (1 + fp')fkfk'] 

+ 2 6{po + Ep, -Uk + u^k') [~fp'{l + fk)fk' - (1 + U')fk{l + fk')] 

+ 8{pq + Ep, +Uk + uJk') Ifp'fkfk' - (1 + /p')(l + /fc)(l + /fc')] 

+ 5(po - Ep> -UJk- UJk') [(1 + /p')(l + /fc)(l + fk') - fp'fkfk'] 

+ 5(po - Ep' +ujk + UJk') [(1 + fp')fkfk' - /p'(l + /fc)(l + /fc')] 

+ 2 Sipo - Ep. +ujk- UJk') [(1 + fp')fk{l + /fc') - fp'{l + /fc)/^] } . (49) 

The function fp was introduced here, 

U = ^E^^ ' (50) 

thermal x particle distribution function being written as fk,fk'- The phase space 

factor here is 

dll„ = - — -; , dUk = -, — T7, . (51) 

^ (27r)3 2Ep' (27r)3 2tUfc ^ ' 

An important point to note is that this spectral function r_ is defined even off the 

mass shell pq ^ \/p^ + M^. Since the negative po region gives r_(— j9o,p, oo) = 

— r_(po,P, C)o), one should be careful in associating indivisual terms of eq.(|49|) with 

physical processes. 

The correlator to 0[A^] is obtained using this r_ in eq. (|45|) along with (^31). This 

result can be compared to a calculation using the thermal field theory in [^. The 

imaginary-time formalism gives the correlator to 0[A^], 

(V3(ri,fi)v3(r2,f2))/3 

A /"^ 4 

= A^{xi-X2)--J d yA^(0)A^{y-xi)A^{y-X2) 

A^ /•z? 
+ ^ / c?V d^y2 A^{yi - y2)A^{yi - y2)A^{yi - y2)A^{yi - Xi)A^{y2 - X2) 
I Jo 

A^ /-/S 

+ xi d^y^d^y^^vi^)^xiyi-y2)^xiyi-y2)\{yi-xi)\{yi-x2) 

+ T i "^'^^ "^'^^ ^x(0)Ax(0)^^(2/i - ^2) A^(yi - xi)A^{y2 - X2) (52) 

^ J-0j^Tj2e~--'^^^-^^^e'P^<'^-'^^F{uj^,p). (53) 

11 



where the variable r is the Euchdean time defined in the range of ~ /5, and uJn = 
27iin/P (n = , ±1 , ±2 ■ ■ ■). The quantities A<^ and A^ are the thermal propagators 
of the (f and x particles, respectively. Each term in eq.(^2|) corresponds to respective 
diagram in Fig. 2. The function F here is calculated as 

F{u:^,p)-' = -ul+p' + M\T) + r duj'-^^^ , (54) 

J—oo UJ — UJ^ 

'' '^''^>"v- r,,tl(^jL^,^^^^] . ,55) 



327r3' 27r ' Jo ul \{1 - e-'^'^^y /3wfc(e/3^fe - 1) 



valid at m <^ T <t^ M where Uk = yk"^ + m'^. The temperature dependent mass 
M'^{T) has correction both of A and A^ orders, but the 0[A^] term is Boltzmann 
suppressed. In the function F the proper self-energy given by a r_ integral is summed 
up to get a geometric form, which is allowed to 0[A^]. The discrete energy sum 
here may be converted by a contour integration to the form given by the real-time 
correlator, which is nothing but eq.(^) along with (^3]). Thus, two methods give 
the same result. 



5 Quantum kinetic equation 

Derivation of the kinetic equation goes in a few steps. We shall only sketch the 
main point of this sequence of arguments. It starts from the self-consistency equation 
for the quantity a{k,t), the Fourier transformed correlator. We first note that the 
odd part of this quantity 

a^{k,t) = a{k,t)-a{-k,t) (56) 

slowly varies. Thus, to order of A^ one can use the infinite time limit; 

a.{k,t) = a-{k,oo) + 0[X^]. (57) 

The basic reason for this simplicity is due to a symmetry of certain terms of eq. 
under Xq ^^ Hq. 

The even part, or a more convenient quantity 

^(*-')-^ '^* 
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which is constrained by 

T{ko,k,t) + T{-kok,t) = -l, (59) 

is used for derivation of the kinetic equation. Under the condition (^) one may use 
the "constant" cr_{k , oo), which can be replaced by 

2nH{k,oo) = a_{k, oo). (60) 

The kinetic equation is derived from the self-consistency relation symbolically 
written as 

r(t)=T[r(t),t], (61) 

where the explicit and implicit (via the function r) time dependence is written. The 
right hand side is understood as a functional of r(t). By dropping higher order OfA"^] 
terms, we may derive from this equation 

dr 

- = -r(r-T[r(t),oo]), (62) 

r = -^ln(T[r(t),t]-r[r(t),oo]). (63) 

The neglected term is dT/dt of order A^ and it was replaced as 

or 

l--^l, (64) 

in the left hand side. One explicitly works out this symbolic equation for our case 
and in the end derives the quantum kinetic equation in the form, 

dr{p,t) _ A^ r (Pki r cPk2 r cPp' f^ dp'p 4 

dt ~ 2po J {2n)^2uJk, J (27r)32cufc, J (27r)3 Jo 2ti ^ ""' ""'^^ ' "^^ 
I 5(4)(p + p'-k,- k,) [rr'(l + A)(l + h) - (1 + r)(l + r')hh] 

+2,^(4) (j9 + p' + k,- k,) [rr7i(l + /2) - (1 + r)(l + r')(l + /O/a] 

+5^'\p + p' + k, + k2) [rr'hh - (1 + r)(l + r')(l + /i)(l + /s)] 

+5(4)(p _ p' _ A;i - k^) [r(l + r')(l + /i)(l + /2) - (1 + r)r7i/2] 

+2<5(^)(p -p' + k,~ k,) [r(l + r')/i(l + /2) - (1 + r)r'(l + f,)f,] 

+5(4) (p _ p' + A;i + A;2) [r(l + r')/i/2 - (1 + r)r'(l + /i)(l + h)] } . (65) 

In the right hand side r = T{p,t),T' = t{p' ,t). In the final step the Markovian 
approximation for the rate was used; the exponential decay law for the quantity 
T[T{t) ,t] was assumed with the pole approximation. 
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The kinetic equation thus derived has a structural resemblance to the conven- 
tional Boltzmann equation. There are however a number of differences. The most im- 
portant one is that there is no mass shell condition for the ip particle; Pq ^ — p ^ 7^ M^. 
Accordingly the function r is a function of the 4-momentum p. Related to this is that 
even the processes not allowed by the energy-momentum conservation on the mass 
shell all contribute to the collision term in the right hand side, for example 1 ^-> 3 
process such as y? ^^ VXX- Even if the equilibrium solution r(p, 00) = l/[e^^° — 1) 
discussed in the preceeding section coincides with the familiar Bose-Einstein form 
at po = Vp^ + ^'^) the pq variable here is defined all the way from —00 to 00. 
These are important differences, although the structural resemblance is impressive 
and suggests deeper understanding. 

Once the kinetic equation is solved, one can compute physical quantities using 
the function r. For instance, the ip energy density is calculated taking the coincident 
limit of two-point correlators, with due consideration of renormalization of composite 
operators. When Fourier transformed, the correlator is given by 



Hence 



a{p,t) = a^{p,t)T{-p,t)=2iTH{p,t)T{-p,t). (66) 



1 1 M^ 

{Ti^) = (-0^ + -{VipY + -^v"^ - (counter terms)) 



1 r d^p 

2 J '(2^ 
S'p 



J 7^^ {pI + Ep}H{p , t)r(p, t) - (counter terms) 
dpo{pl + EnH{p,t)T{p,t) 



(27r)3 Jo 

1 f d^T) f°° 
+ 2 y 72^ J ^^° ^^0 + Ep}H{p , t) - (counter terms) . (67) 



(2vr) 
This is further simplified using H{p , t) ~ H{p , 00) and 

1 f d^p f°° j_ r^{p, 00) , 1 , ^r.3n 



^po T Z.4 + + QM- (6^ 



2 J {2Txfk "^"{po + E.y 2 



Thus, 



J^l dpo{pl + El}H{p,oo)r{p,t) + c\' — 
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^ 1.4 X 10-^ (70) 



69120 

dropping Boltzmann suppressed 0[A^] terms. In the last step we replaced the spectral 
weight H by the delta function, since perturbative analysis works for this quantity. 
In this computation we assumed the mass hierarchy for the light x particle, m ^T. 
We used in this derivation the formula (E9|) to get {T <^ M assumed) 



1 f d^p /""^ r_(p,oo) 1 



2X^ J dUk J dUk> {P + k' 2) /,/,, j dUp J rfn, 



{2Trf5''^\p + p' + k + k') 



r-^ 



^-^2] dll>.J dU,,{e + k'^)f,f,,l dp^^=cX'—, (71) 

dropping Botzmann suppressed minor terms. Physical processes that contribute to 
this equilibrium result are 

XX^^^' , X^ ^^'x , (72) 

(p' being the partner ip particle. All other channels are suppressed by the Boltzmann 
factor. 

The time dependence of physical quantities appears via the function T{Ep,p,t) 
in eg . (|69|) . This function defined on the mass shell {Ep = y/p'^ + M^) follows the 
usual Boltzmann equation, as seen from the kinetic equation (^) along with the 

approximation, 

2n 
a-{p,oo)^^ 6{po - Ep) . (73) 

It is thus clear that in a state near thermal equilibrium T{Ep,p,t) ~ l/{e^p^'^ — 1). 
Hence for T/M < l/y^\ ^ 3 x lOVA 

The same equilibrium result was derived using the thermal field theory in [^ . Deriva- 
tion given here is a separate and independent confirmation of the equilibrium density 
from the kinetic approach. The temperature power thus derived (oc T^) differs from 
that given in ref. 0, which should be corrected. 
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As discussed in [0, the time evolution equation in cosmology is approximately 
given by (at m <^ T <^ M) 

^ + 3i7n, = -M«-ny , (75) 

MT T^ 

,,^^..(MLf,^,-MIT^^^.T_^ (76) 



where H is the Hubble rate ~ LlyNT"^ /m^i. Numerical analysis shows that the 
picture of the sudden freeze-out holds such that until the freeze-out the num- 
ber density follows the equilibrium value rieq. Thus, the most important part of 
cosmological application is to estimate the freeze-out temperature using the new 
equilibrium number density. This estimate is done by equating the annihilation rate 
{av) ricq to the Hubble rate H. 

The close connection with the equilibrium thermal field theory is made evident 
by taking the infinite central time limit; the non-equilibirum correlator defined here 
approaches the equilibrium value given by the thermal field theory, 

{ip{x)ip{y)) = tr Pi ip{x)ip{y) -^ tr e"^^'°* ip{x)ip{y)/tT e"^^'"* . (77) 

The equilibrium correlator is time-translation invariant, hence the central time de- 
pendence disappears for the correlator and the non-equilibrium correlator recovers 
this invariance. Moreover, our formula for the correlator and physical quantities such 
as the energy density indicates that the approach towards thermal equilibrium de- 
scribed by the first term in eq.(|69D is essentially determined by a thermally averaged 
rate of on-shell quantities, consistent with the previous estimate of this quantity in 
derivation of the kinetic equation. 

Thus, the most important change for cosmological estimate of the relic mass 
density is the new equilibrium quantity. This equilibrium value is precisely given 
by the Gibbs formula, using the total hamiltonian, not its free part. This is the 
main reason the use of thermal field theory in our companion paper 0] is justified in 
computation of the relic mass density. It has to be kept in mind that the simple ideal 
gas form of distribution function must be questioned in the low temperature region. 
Even a higher order term in coupling is more important than the zero-th order 
Boltzmann suppressed term when the Boltzmann suppression factor is huge. The 
appearance of our temperature power term for the equilibrium number density is due 
to a continuous integral of these Boltzmann suppressed terms which dominates over 
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a single contribution of quasi-particle object given by approximating a Breit-Wigner 
form of integral by a single complex pole. 

Despite this rather obvious change, the thermally averaged Boltzmann equation 
is hard to justify at low and intermediate temperatures, or its slight modification 
is difficult to lead to the relevant correction. Although our work justifies in some 
sense the Boltzmann equation at high temperatures, it should be kept in mind that 
full quantum mechanical treatment is essential to derive the correct low temperature 
result. 
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Figure caption 

Fig.l 

Diagram for the spectral weight of two hght x particle intermediate states denoted 
by broken lines. 

Fig.2 

Diagram for thermal correlator to 0[A^]. The solid and broken lines correspond 
to the if and x particle propagators in thermal equilibrium, respectively. 
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Fig. 1 
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